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I summarise what recent lattice calculations tell us about the large-A'^ limit of 
SU(Af) gauge theories in 3+1 dimensions. The focus is on confinement, how close 
SU{oo) is to SU(3), new stable strings at larger N, deconfinement, topology and 
S-vacua. I discuss the efTective string theory description, as well as master fields, 
space-time reduction and non-analyticity. 



1. Introduction, Overview, Conjectures 



It is 30 years since it was proposed that it might be useful to think of QCD 
r> I as a perturbation in 1 /N around the N = oo theory ^ . As is apparent from 

' the talks at this meeting, this has been a very fruitful idea ^. However we 

still do not have a quantitative control of the SU(A^ — oo) theory and the 
phenomenology needs to assume, for example, that it is confining and, of 
course, 'close to' SU(3). Lattice simulations can attempt to answer such 
questions directly (albeit never exactly) and there has been substantial 
progress in doing so this last decade, first in D=2+l dimensions ^, which I 
do not discuss here, and then in the physically interesting case of D=3+l. 
(For a review of work in the earlier 80's, when lattice calculations were 
not yet precise enough to be so useful, see e.g. **.) Here I focus on what 
these 'modern' lattice calculations teach us about the properties of SU(A^) 
gauge theories at large N. I will begin with some motivation for these 
calculations. 

A gluon loop on a gluon propagator comes with a factor of g'^N. One 
easily sees that is in fact the smallest power of the coupling that comes 
with a factor of N. So if one wants an A'^ — > oo limit that is not given by 
either a free field theory or by infinite order diagrams on all length scales (in 
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neither case would we get something hke QCD), then one needs to take the 
hmit keeping g^N fixed ^. In Section 6 we shall see that this standard all- 
order perturbative statement also appears to hold non-perturbatively. Us- 
ing 't Hooft's double-line notation for gluons, diagrams can be categorised 
as lying on surfaces of different topology, with more handles corresponding 
to higher powers of 1/-/V, so that in the N = oo limit only planar diagrams 
survive ^. As the coupling g^N becomes strong, the vertices of the dia- 
gram fill the surface more densely, defining a world-sheet of the kind one 
might expect in a string theory. This suggests that perhaps as A?^ ^ cxi 
the gauge theory can be described as a weakly interacting string theory ^. 
If the theory is, as one expects, linearly confining - and in Section 3 we 
present evidence that it is - then the confining flux tube will behave like 
a string at long distances, described by some effective string theory. At 
large A'' this will presumably coincide with the string theory that describes 
the SU(A'') gauge theory. In Section 3 we shall discuss what a study of 
the string spectrum teaches us about this string theory. These 'old' string 
theory ideas ^ have been recently complemented by the realisation that at 
g'^N ^ oo and N ^ oo gauge theories have a dual string description that 
is analytically tractable ^. Determining the effective string theory numer- 
ically should provide useful hints about what the dual theory might be in 
the physical weak-coupling limit, g-^N ^ 0. 

In a confining theory there are no decays at A'' = oo. This is in contrast 
to what would happen in a non-confining theory where a coloured state 
could have a finite decay width into other coloured states. But once we 
constrain states to be colour singlet we reduce the density of final states 
by factors of N so that all decay widths vanish. In addition there is no 
scattering between the colour singlet states. Think of two propagating 
mesons. A meson propagator is like a closed quark loop. Exchange two 
gluons between these two closed loops and you clearly gain (up to) a factor 
of N, but at the cost of g** oc So, no scattering at TV = oo. However 

it is also easy to see that within a single closed loop, planar interactions 
give factors of g^N, are not suppressed, and so there are non-trivial bound 
states. So we have what looks like a free theory, but it has a complex bound 
state spectrum and so is non-trivial. If one is going to find room in D=3+l 
for notions of e.g. integrability, it is here in the N = oo limit that they 
might find a suitable home. 

Is this theory with no decays and no scattering similar to the observed 
world of the strong interactions? First we need the SU(oo) theory to be 
linearly confining - and we provide evidence for this in Section 3. Now we 
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can ask: is this confining theory close to SU(3)? In Section 4 we calculate 
the lightest masses in the spectrum for several values of N, and find that the 
SU(3) mass spectrum is indeed very close to the (extrapolated) spectrum 
at A'' = DO. This provides support for the phenomenological relevance of 
the SU(oo) theory. And this provides motivation for trying to understand 
that theory much better. In Section 3 we also see that the effective string 
theory describing long confining flux tubes appears to be in the bosonic 
string universlity class. More surprisingly, the energy of shorter strings is 
close to the Nambu-Goto prediction and at smaller N, where this question 
can be addressed, the string condensation temperature is very close to that 
of the Nambu-Goto string action. A new phenomenon for AT > 3 is the 
existence of new stable strings. In Section 7 we summarise the latest lattice 
calculations of the corresponding string tensions and find they lie between 
the 'Casimir scaling' and 'MQGD' ^ conjectures. We remark how these 
new strings can contribute to an A/'-dependence of the mass spectrum even 
in the confining phase, contrary to naive expectations. In Section 8 we learn 
how the large- gauge theory deconfines. Contrary to some speculations 
that the N = oo transition might be second order, partly motivated by 
the weakness of the first order transition in SU(3), we will show it to be 
robustly first order. 

At AT = oo the expectation value of a product of gauge invariant op- 
erators factorises into the product of the respective expectation values - 
by the same argument that there is no scattering. This suggests that a 
single gauge orbit - Witten's Master Field ^ - dominates the Path Inte- 
gral calculation of all the physics in the confined phase. Since the physics 
is translation invariant, so must the Master Field be (for gauge invariant 
quantities). This suggests that all we need to know is the field in an arbi- 
trarily small region to know it everywhere - even on one point if that can 
be made precise by a suitable regularisation. On the lattice this is achieved 
through (twisted) Eguchi-Kawai reduction 

In the deconfinement transition one sees explicitly how the large- A^ be- 
haviour of various quantities - latent heat, interface tension, fluctuations - 
means that a 'phase transition' occurs on ever smaller volumes as A^ ^ oo. 
We shall also see, by a heuristic but physical argument, why the imposition 
of twisted boundary conditions is required to remain in the right phase. 
Precisely at the deconflning temperature, Tc, there is a different Master 
Field of the Euclidean Path Integral for the confined and deconfincd phases. 
Through hysteresis this extends either side of this temperature; conceivably 
to all T. Indeed the Euclidean system possesses N different phases, and 
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hence Master Fields in the deconfincd phase. This multiphcation of master 
fields is not peculiar to deconfinement. For example, intertwining 0- vacua 
would lead to N non-degenerate vacua at ^ = which become abso- 
lutely stable at A'' = oo, as discussed in Section 9. Each of these vacua will 
have its corresponding master field. 

For > 5 one finds, in the lattice gauge theory with the standard 
Wilson plaquette action, a first order 'bulk' phase transition at a particular 
value of the inverse 't Hooft coupling Ac(a) = g'^{a)N ~ 1/3. (We write 
the bare coupling as a running coupling on the scale a.) This is essentially 
the same as the £' = l + l,iV = oo Gross- Witten phase transition . 
For A(o) > Ac(a) the vacuum is non-perturbative on all lengths scales, 
so that the confining string tension is 0(1) in lattice imits. For A(a) < 
Ac (a) the vacuum is perturbative on the shortest distance scales, and hence 
asymptotically free as o — > 0, and the string tension is 0(1) in physical 
units. This transition appears as a lattice peculiarity but, as has recently 
been discovered I4,i5^ there appears to be an analogous non-analyticity at 
N = oo that occurs as we increase the size of a Wilson loop: at a certain 
critical size the eigenvalue spectrum of the loop changes non-analytically 
14,15^ This is possible because at = oo the number of physically relevant 
degrees of freedom per unit volume is infinite. It appears that this critical 
size is fixed in physical units and will survive in the continuum limit. 

There are, of course, other non-analyticities as iV — > oo. For example, 
we shall see in Section 9 that the instanton size distribution exactly van- 
ishes for sizes up to some critical size. However this can be understood as 
due to the factor exp(— Stt^/^'^) oc exp(— cA'') that dominates the weighting 
of small instantons. We do not know of any such simple argument in the 
case of Wilson loops. Indeed we might conjecture that this non-analyticity 
provides an explanation for the puzzlingly rapid transition between short 
and long distance physics that is observed experimentally i.e. as soon 
as one is at values of where one can apply perturbation theory, one 
finds that there is little room for the higher twist operators that one might 
expect to parametrise the transition from perturbative to non-perturbative 
physics - 'precocious scaling'. The non-analyticity discovered in sug- 
gests that for I < Ic we can calculate the wilson loop perturbativcly, while 
for I > Ic it is confining and non-perturbative. At A^ = oo the transition is 
infinitely sharp; in SU(3) and QCD it might become a very rapid cross-over, 
explaining the phenomenon of precocious scaling. 

The presence of such a 'phase transition' as we increase the distance, 
might effectively disconnect the confining theory from its short distance 
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perturbative framework. In this disembodied confining theory the coupUng 
is never small, it is confining on all available length scales and one never 
needs to discuss gluons. This raises, for example, the possibility of a dual 
string theory in which the coupling need not be large. Such a dual theory 
might be analytically tractable. It also raises the possibility that at = oo 
the same confining theory can have different ultraviolet completions. That 
is to say, to solve the theory we do not necessarily need to solve the full 
non-Abelian gauge theory. These conjectures are speculative, of course, 
but they certainly provide motivation for clarifying ^"^ the nature of this 
remarkable non-analyticity. 

2. Lattice 

We will calculate Euclidean Feynman Path Integrals numerically. This 
requires a finite number of degrees of freedom, so we discretise continuous 
space-time and make the volume finite by going to a hypercubic lattice on 
a 4-torus. Since the theory is renormalisible and has a mass gap, the errors 
induced by this should rapidly disappear as the lattice spacing is reduced 
and the volume enlarged. The lattice spacing is a and the size of the fi-toms 
is lij, = aL^. The degrees of freedom are SU(Af) matrices, Ui, defined on 
the links I of the lattice. The partition function is 

Z{(3) = j I]dC/,e-'^^.{i-^^''^»''> ; /3 = ^ (1) 
J ^ 9 

where Up is the ordered product of matrices around the boundary of the 
elementary square (plaquette) labelled by p and is the bare coupling. 
This is the standard Wilson plaquette action and one can easily sec that 
for smooth fields it reduce to the usual continuum gauge theory. Since 
the theory is asymptotically free and since the bare coupling is a running 
coupling on length scale a, the continuum limit is approached by tuning 
[3 = 2N/g'^{a) oo. As we remarked earlier, one expects from the di- 
agrammatic analysis that for large N the value of a is fixed in physical 
units (e.g. in units of the mass gap) if one keeps the 't Hooft coupling 
A(a) = g'^{a)N fixed i.e. f3 cx N"^. This will be confirmed below. 

The lattice path integral in eqn(l) is no easier to calculate analytically 
than the original continuum version. However because the number of inte- 
grations is now finite, we can attempt a numerical evaluation. The number 
of integrations is large and so the natural method to use is the (Markovian) 
Monte Carlo. The Monte Carlo generates 'points' in the integration space. 



Februciry 1, 2008 19:47 Proceedings Trim Size: 9in x 6in 



draft2 



6 

Each such 'point' is an explicit lattice gauge field i.e. an SU(A'') matrix on 
every link of the lattice. These fields are generated with the measure 

W = []dJ7ze-'^^.^i-^^^^"-> (2) 

so if we generate ric such 'points', i.e. {U^{n);ii = 0, ..,3;n = 
1,...,Z/^}^ ; / = l,...,nc, then the expectation value of will be just 
the average over these fields: 

1 1 
(*L(C/)) = -^*i([/^)±0( — ). (3) 

I have made explicit here the statistical error which decreases as the square 
root of the number of field configurations - as one would expect for such a 
probabilistic estimate. 

We calculate masses from Euclidean correlation functions 

C{t = ant) = {cpHmO)) = Kn|</.(0)|mc)|2e-«^""* ; < Ei+, 

n 

(4) 

which we evaluate numerically as just described. Note that all energies will 

be obtained in lattice units, ai?„. At large rit the lightest state will dominate 
C{t = ant) and can be easily extracted. Unfortunately the statistical error 
in the calculation of eqn(3), is more-or-less independent of nt, since the 
average fluctuation squared around the correlator is itself a higher-order 
correlator which, one can easily verify, has a disconnected piece. Thus the 
error to signal ratio grows exponentially with nt and one needs C{t = ant) 
to be dominated by the lightest state at small nt i.e. one needs to be a 
good wave-functional for the desired state. Standard techniques now exist 
to achieve this, and can be used within a variational calculation, based on 
the exp(— aiJ) implicit in C{t = ant), to obtain excited as well as ground 
state energies However it should be apparent that the larger the energy, 
the less accurate the calculation. 

To calculate glueball masses we use operators that are based on con- 
tractible Wilson loops. We calculate the string tension from the energy of 
the lightest flux loop that winds around a spatial torus, and use operators 
based on a Wilson line that encircles the torus. In Fig. 1 1 show an example 
of the latter The calculation is in SU(6) on a 12 x 14^ lattice with the 
flux loop winding around the x-torus. Shown also is the best single expo- 
nential fit (actually a cosh because of the periodicity in t). It is clear that 
it dominates the correlator from very small t - indeed the overlap of the 
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Figure 1. Correlation function where the hghtcst state is a flux loop that winds around 
the spatial torus. Best single energy cosh fit shown. 

operator on the flux loop is ~ 0.97. This is achieved by iterative smearing 
of the fields and by a variational calculation (see e.g. for details). 

3. Confinement and Strings 

Consider one spatial torus of size I and all the other tori large. We calculate 
the mass of the lightest flux loop that winds once around this torus. We 
expect its energy to be 

,,.l^oo J 7r(rf-2)c 
m{l) = al y (5) 

where c (times the dimensional factor d — 2) is the central charge of the 
effective string theory that describes the long-distance properties of the 

confining fiux tube. In Fig. 2 I show the results of a calculation of am{l) 
in SU(6) at a fixed value of a We see linear confinement, and a fit with 
eqn(5) to I > 10,12 gives c = 1.16(8), 1.09(10) respectively. This tells us 
that the effective string theory at long distances is a simple bosonic theory 
in the universality class of the Nambu-Goto action. Indeed, if we fix the 
coefficient of the l/l term to the bosonic string value 7r/3, then we find that 
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Figure 2. The masses of the lightest, and first excited, o, fc = 1 flux loops that wind 
around a spatial torus of length I in the SU(6) calculation at /3 = 25.05. The dotted 
lines are the predictions of the Nambu-Goto string action, as in eqn(7). The dynamical 
lower bound on the string length is Imin = l/iTc — 6.63. 

we can obtain an acceptable fit to our whole range of I: 

TT 19 4f3 2) 

ai;o(/) = 0.06358(20); : l>7 ,xVndf = 0.9 (6) 

This is remarkable: since there is a minimum length for a periodic flux 
loop, which is Imin = l/a^c — 6.63 in the present calculation, the fit in 
eqn(6) essentially works all the way down to the shortest possible strings. 
(If wc fit the 0{l/l) term as well, then its coefficient comes to 0.94(16).) 
Since the corrections in the pure gauge theory to iV = oo are 0{1/N^), we 
can assume that all this is also true of the SU(iV = oo) theory. In physical 
units the lattice spacing is small, a^/a ~ 0.25, so we can assume that this is 
also true of the continuum limit. Finally, the longest string is 20a ~ S/y^cr 
so we can assume we are seeing the asymptotic behaviour of a long string. 
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In the Nambu-Goto string theory the spectrum is given by 



21,22 



d-2 



)}' 



(7) 



n — 



24 



In Fig. 2 we show that the one parameter fit with Eq{1) to the hghtcst 
string mass works not too badly all the way down to the minimum possible 
string length, 1^ = 1/Tc- The is too large to be acceptable, but is much 
smaller than for SU(4). This leaves open the intriguing (and unexpected) 
possibility that the Nambu-Goto string action describes confining strings 
on all length scales at = oo. 

4. Spectrum 

The lightest J^*^ = 0++ and 2++ gluebahs turn out to be the lightest 

states in the a — 0,V = oo SU(A'") gauge theory. In Fig. 3 I plot the ratios 
of these masses to the (simultaneously calculated) string tension, obtained 
after a continuum extrapolation of the lattice results for each value of N 
We see that a modest 0{1/N'^) correction suffices to fit the ratios for 
N > 3: for these quantities SU(3) is indeed close to SU(oo). 



This is for only two states of course. For a much more detailed compar- 
ison of SU(3) and SU(8) see ^3. 



5. Pomeron 

In ^3 you will also find estimates for higher spin glueballs. This requires 

novel lattice techniques because of the reduced cubic rotational symmetry 
of our lattice. This calculation enables us to ask if the Pomeron trajectory 
is the leading glucball trajectory (ignoring mixing). Fig. 4 contains a Chew- 
Frautschi plot for the PC = ++ sector of the D=3-|-l SU(3) gauge theory. 
Recall that l/2TTa ~ IGeV is roughly the slope of the usual mesonic Regge 
trajectories. Assuming linear trajectories, the best fit to the intercept ao 
and slope a' of the leading trajectory is 



This provides quite convincing evidence for the old conjecture that the 
Pomeron is in fact the leading glueball trajectory. It is interesting to note 
that in 2+1 dimensions the leading glueball trajectory has a very low in- 
terecept and so is unimportant at high energies 23.24_ -jjiig jg -(^iig (.^gg uQ-j^ 
only for SU(3) but for larger N as well 



27r(Ta' = 0.281(22) 



ao = 0.93(24) 



(8) 
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Figure 3. The lightest 0+^", •, and 2++, o glueball masses expressed in units of the 
fundamental string tension, in the continuum limit, plotted against l/N'^. Dotted lines 
axe extrapolations to N = oo. 




Figure 4. Chew-Frautschi plot of PC = ++ states in the continuum SU(3) gauge 
theory. The leading Regge trajectory is shown. 
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6. Coupling 

The diagrammatic expectation is tliat one keeps fixed the 't Hooft coupling 
A = g^N for a smooth large-A^ limit. For theories with running couplings 
(as here) this becomes 

Xil^) ^ {I ^=°° hidoi N (9) 

where the length scale / is expressed in units of ^, which is some physical 
mass in the SU(7V) gauge theory that one expects to have a smooth, non- 
zero large- -/V limit. We expect that if we plot A(Z/x) against then as 
— > oo all the curves for different N will lie on top of each other. 
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Figure 5. The square root of the string tension in lattice units, a^cr, plotted against 
the 't Hooft coupling, \j = g^N. 

Let US choose fi — y/a and / = a. Then we can use (3 = 2N / to define a 
coupling g^{a). This is not ideal because of large lattice spacing corrections. 
It is well known that the 'mean field improved' coupling, (3{ReTTup)/N — 
2N/g'j{a), is much closer to standard definitions. Defining A/ (a) — g'j{a)N 
we can plot its value against the value of a^/a obtained in the calculation 
with the corrsponding (3. We plot the results in Fig. 5 for 2 < < 5. 
This provides nice non-perturbative confirmation of the usual perturbative 
result. 
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7. k-strings 

Sources that transform a.s ip ^ z^ij) under a gauge transformation of the 
centre, z € Zjv C SU{N), cannot be screened by gluons to a source that 
transforms otherwise. Thus each sector k has its own stable string tension, 
(T/j. Lattice calculations have shown these are bound, ah < ka, and 
have focusscd on a range of values that includes the Casimir scaling and 
'MQCD' conjectures ^. Results from the most recent work on this topic 
are listed in Table 1. The results lie between the Casimir Scaling and 
MQCD predictions. 





CTfc/ 


a 




(N,k) 


Casimir scaling 


this paper 


'MQCD' 


(4,2) 


1.333 


1.370(20) 


1.414 


(4,2) 


1.333 


1.358(33) 


1.414 


(6,2) 


1.600 


1.675(31) 


1.732 


(6,3) 


1.800 


1.886(61) 


2.000 


(8,2) 


1.714 


1.779(51) 


1.848 


(8,3) 


2.143 


2.38(10) 


2.414 


(8,4) 


2.286 


2.69(17) 


2.613 



The values in Table 1 are obtained after an extrapolation to the contin- 
uum limit: 

^(a) = ^(0)+caV (10) 
a (7 

from a range of a where the leading correction proves sufficient. (It is 
known that the leading correction for the plaquette action is 0(a^).) I 
show in Fig. 6 an example from SU(4), which superimposes two sepa- 
rate calculations. It is clear that the continuum extrapolation is very well 
determined. This is also true of the glucball calculations discussed earlier. 

All this has interesting implications for the mass spectrum. One expects 
that (most) highly excited glueballs can be described as closed loops of flux 
i.e. closed strings. Normally one thinks of flux tubes in the fimdamental 
representation. However closed loops of A;-strings will do just as well. These 
will scale with ak- Thus as N increases we expect the mass spectrum to 
acquire new towers of states, based the stable strings with k < N/2. These 
towers will be copies of each other scaled up from the fundamental by 
^/ok/a. There is a hint of such a state in the comparison between SU(3) 
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0.125 0.15 



Figure 6. Ratio of the tension of the k = 2 string to tliat of the fundamental (fc = 1) 
string in SU(4), for the anisotropic (o) and isotropic (•) calculations. Plotted against 
the square of the (spatial) lattice spacing. Continuum extrapolations are shown. 

and SU(8) in . If follows Casimir scaling then this will lead to an 
'unexpected' 0(1 /N) variation in some mass ratios. It also implies that 
the 'entropy' in the confined phase is not TV-indcpcndcnt, in contradiction 
to usual assumptions (although this will not upset usual conclusions). 



8. Deconfinement and String Condensation 

It is well-known that the SU(2) deconfining transition is second order and 
that SU(3) is 'weakly' first order. Recent numerical studies of deconfine- 
ment in SU(Af > 3) gauge theories have shown that the transition be- 
comes more strongly first order as iV t and it is clear that it is robustly 
first order in the N = oo limit. 

In Fig. 7 I plot the value of the deconfining temperature in units of the 
string tension, for 2 < < 8. We see that a modest 0{1/N^) correction 
suffices to fit these values. This is perhaps puzzling given that SU(2) is 
second order. 

If nothing else happens as one increases T, one expects a string con- 
densation deconfining transition to occur simply because for strings of 
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Figure 7. The deconfining temperature in units of the string tension for various SU(JV) 
gauge theories. Large N extrapolation shown. 

length I the thermodynamic suppression exp{—E/T} ^ cxp{—al/T} is out- 
weighed by the number of strings oc expl+Cg/} once T > a/c = T^, where 
we call the string condensation temperature. Thus any second order 
deconfining transition may be driven by string condensation and the value 
of T^/a will then tell us something about the effective string theory. For 
SU(iV) gauge theories we have 

SU{2),d = A 

SU{2),d=3 (11) 
SU{3),d=3 

and it is interesting to compare it to the Nambu-Goto value 

11= I ^ f°-^^^ ■ ^ = '^ (12) 

Y(d-2)7r |o.977 : d = 3 ^ ' 

obtained from the vanishing of Eq{1) in eqn(7). 

As expected, one finds that the latent heat Lh oc and that the 
interface tension, (Ted, grows with TV. It is easy to see that if acd « TV 
then at = oo the transition becomes infinitely sharp, i.e. a real phase 
transition, even on a finite volume and there is no hysteresis. If N'^ 




0.709(4) 

1.12(1) 

0.98(2) 
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then the hysteresis may be large - perhaps even infinite! In that CclSG ^ clS 
we increase T past Tc in the confined phase, we may continue to remain 
confined to a high enough T that we encounter the string condensation 
transition T*. If so, then we can study string condensation in the interesting 
N = oo Hmit. 

Although one finds that one can calculate Tc on ever smaller volumes 
as A'' I, one also knows that once one of the spatial sizes drops below 
^ 1/Tc, flux loops that wind around that torus will condense, and this 
spontaneous symmetry breaking will alter the thermal physics. This looks 
like an obstacle to reducing the volume further as iV "f . However one knows 
how to prevent this - twisted boundary conditions can force a domain wall 
that restores the symmetry For small sizes the wall unravels and the 
system sits in the symmetric core. Then one can continue to reduce V with 
A''. This provides a physical interpretation of twisted Eguchi-Kawai. 

9. Topology, Instantons, 0-vacua 

The topological susceptibility Xt has a non-zero limit at A = oo with a 
value that fits in with the Wittcn-Vcncziano analysis of the rj' mass 25, 30^ 
This is of course a non-leading effect -''^ since (Q^) — {Q)^ = to leading 
order. 

In lattice calculations a very small instanton is a very big spike in the 
FF density and is easy to see. So in the instanton density plot in Fig. 8 
the rapid suppression of small instantons with A' is reliable - and of course 
expected from the 



factor that dominates the weight as p ^ 0. This has consequences for the 
numerical calculation. Since the Monte Carlo deforms the fields 'continu- 
ously' (one link matrix at a time) the only way to change Q is to shrink an 
(anti)instanton out of the lattice (or the reverse). Thus the probability to 
change Q is linked to the probability to find an instanton of size p ~ fewa 
and so is suppressed as oc a^^^/^~^. So for large A' a lattice calculation 
gets stuck in a fixed topological sector "^^ . You could then try to calculate 
{Q^) from a large sub-volume of a very large lattice volume. 

Larger instantons are more ambiguous but Fig. 8 certainly suggests the 
interesting possibility 
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Figure 8. The instanton size density, -D(p), for N = 2(*), 3(+), 4(o), 8(») on 16"* lattices 
with a ~ l/8Tc. 



Usual counting tells us that the physics should be a function of 6/N 
and so periodic in 2-jtN. However Q integer tells us it is periodic in 2it for 
all TV. This implies that there arc in fact vacua that intertwine as 9 
varies. At 9 = there are N 'vacua' which become stable (exponentially 
fast) as A/' — > 00. One can see these as a continuation from J\f = 1 SUSY of 
the N degenerate vacua associated with gluino condensation As a first 
step, on the lattice side, it has been shown that the moments of Q at 
^ = do suggest a vacuum periodicity of 2TrN rather than 2tt. What one 
would like to see are the near-stable 'vacua' directly. It is easy to calculate 
their energies 

En{e = 0) = lxt{'2nnfV ; AT ^ oo (15) 

where Xt — 0.39-^(7 ~ 2Tc/3. Thus sometimes one might tumble into an 
n^O vacuum from the dcconfined phase at Tc, but it should be easier from 
the bulk phase. This could be done on the lattice. 
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10. Some directions ... 

It would be very nice to have the spectrum of QCDm=oo with light quarks: 
glueballs and mesons, ground state and excited, all stable and non-mixing. 
Since this can be approached through extrapolating quenched theories in 
N, this is feasible. 

As I pointed out, the intertwining 0- vacua might well be directly acces- 
sible in lattice calculations. It would be nice to follow them back to = 1 
SUSY by varying, say, a gluino mass parameter. 

Space-time reduction and the non-analyticities at A'' = oo might have 
important phenomenological and theoretical implications as conjectured 
earlier in this talk. 
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